Optical Sum Rule in Strongly Correlated Systems 
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We discuss the problem of a possible "violation" of the optical sum rule in the normal (non super- 
conducting) state of strongly correlated electronic systems, using our recently proposed DMFT+E 
approach, applied to two typical models: the "hot - spot" model of the pseudogap state and dis- 
ordered Anderson - Hubbard model. We explicitly demonstrate that the general Kubo single band 
sum rule is satisfied for both models. However, the optical integral itself is in general dependent 
on temperature and characteristic parameters, such as pseudogap width, correlation strength and 
disorder scattering, leading to effective "violation" of the optical sum rule, which may be observed 
in the experiments. 

PACS numbers: 74.25.Gz, 71.10.Fd, 71.10.Hf, 71.27.+a, 71.30,+h, 74.72.-h 



I. INTRODUCTION 

Many years ago Kubo [lj has proven the general sum 
rule for diagonal dynamic (frequency dependent) conduc- 
tivity (t(u>), which holds for any system of charged parti- 
cles irrespective of interactions, temperature or statistics. 
This sum rule is usually written as: 
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Rea(uj)duj = 



(1) 



and 



where r specifies the type of charged particles, n r 
are the respective densities and charges. 

For the system of electrons in a solid Eq. {!]) takes the 
form: 



Rea{ui)dui = —j- 



(2) 



where n is the density of electrons and uj pl = 47r ™ is the 
plasma frequency. 

However, in any real experiment we are not dealing 
with an infinite range of frequencies. If one considers 
electrons in a crystal and limits himself to the electrons in 
a particular (e.g. conduction) band, neglecting interband 
transitions, the general sum rule reduces to the single 
band sum rule of Kubo [l], Q : 



W = 



Rea(uj)dto = f(uj c 
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(3) 



where e p is the bare dispersion as defined by the effec- 
tive single band Hamiltonian, while n p is the momen- 
tum distribution function (occupation number), which is 
in general defined by the interacting retarded electronic 
Green's function G R (e,p) [3], [1]: 



den(e)ImG R (e, p) 



(4) 



where n(e) is the usual Fermi distribution. In Eq. © 
u) c represents an ultraviolet cut-off, frequency, which is 
assumed to be larger than the bandwidth of the low en- 
ergy band, but smaller than the gap to other bands. The 
function f(to c ) accounts for the cut-off dependence, which 
arises from the presence of Drude spectral weight beyond 
u) c Q and is unity if we formally set ui c to infinity while 
ignoring the interband transitions. 

Although the general sum rule is certainly preserved, 
the optical integral W(ui c ,T) is not a conserved quan- 
tity since both f(ui c ) [f| and n p 0, [1] depend on tem- 
perature T, and also on details of interactions This 
dependence of W on T and other parameters of the sys- 
tem under study has been termed the "sum rule viola- 
tion" . It was actively studied experimentally, especially 
in cuprates, where pronounced anomalies were observed 
both in c-axis and in-plane conductivity both in normal 
and superconducting states [1, H, [l(| El, 0) EH • 

The finite cut-off effects were extensively studied in 
several theoretical papers on the T dependence of the 
optical integral 0, H, 0|- In Refs. [1, 0], the effect 
of the cut-off was considered in the context of elec- 
trons coupled to phonons. In a simple Drude model, 
a(uj) = (Wpj/47r)/(l/T — iui) and the sum rule can only 
be "violated" due to the presence of f(u> c ). Integrating 
over tu and expanding for lo c t >> 1, one can see that 
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For infinite cut-off, /(w c ) = 1 and W = Wp//8, but for 
a finite cut-off f(u> c ) contains the term proportional to 
1/uj c t. If 1/r changes with T, then one obtains a sum 
rule "violation" even if to p i is T independent [a, [7| . Other 
aspects of cut-off dependence were discussed recently in 
detail in Ref. @. 

In the present study we neglect the cut-off effects in 
optical integral from the outset. Our goal is to study 
W dependence on T and a number of interaction param- 
eters, determining the electronic properties of strongly 
correlated systems, such as cuprates. In this context we 
shall discuss the problem of a possible "violation" of the 
optical sum rule in the normal (non superconducting) 
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state of strongly correlated electronic systems, using our 
recently proposed DMFT+E approach [3, El, as 
applied to dynamic conductivity in two typical models of 
such systems: the "hot - spot" model of the pseud ogap 
state [19( and disordered Anderson - Hubbard model 20] . 
Our aim is both to check the consistency of DMFT+E ap- 
proach as applied to calculations of optical conductivity, 
as well as to demonstrate rather important dependences 
of the optical integral W not only on T, but also on such 
important characteristics as pseudogap width, disorder 
and correlation strength, making (single band) sum rule 
"violation" rather ubiquitous for any strongly correlated 
system, even neglecting the cut-off effects. 



II. OPTICAL SUM RULE IN THE 
GENERALIZED DMFT+E APPROACH 

Characteristic feature of the general sum rule as ex- 
pressed by Eqs. ([3]), ^ is that the integral W over fre- 
quency in the l.h.s. is calculated via two-particle prop- 
erty (dynamic conductivity, determined by two-particle 
Green's function, in general, with appropriate vertex cor- 
rections), while the r.h.s. is determined by the single- 
particle characteristics, such as bare dispersion and oc- 
cupation number ((4|) (determined by a single-particle 
Green's function). Thus, checking the validity of this 
sum rule, we are in fact thoroughly checking the con- 
sistency of any theoretical approach, used in our model 
calculations. 

Our generalized dynamical mean field theory 
(DMFT+S) approach [3 E! El, supplying the 
standard dynamical mean field theory (DMFT) [I?], El 
with an additional "external" self-energy E (due to any 
kind of interaction outside the scope of DMFT, which is 
exact only in infinite dimensions), provides an effective 
method to calculate both single- and two-particle 
properties [3 HI. The consistency check of this new 
approach is obviously of great interest by itself. We 
shall also see, that it gives a kind of a new insight in the 
problem of sum-rule "violation" . 



A. Pseudogap state, the "hot spots" model 

Pseudogap phenomena in strongly correlated systems 
have essential spatial length scale dependence [21(. To 
merge pseudogap physics and strong electron correla- 
tions we have generalized the dynamical-mean field the- 
ory E3, El| by inclusion of the dependence on correla- 
tion length of pseudogap fluctuations via additional (mo- 
mentum dependent) self-energy E p (e). This self-energy 
E p (e) describes non-local dynamical correlations induced 
either by short-ranged collective SDW-like antiferromag- 
netic spin or CDW-like charge fluctuations [13, HI . 

To calculate E p (e) in two-dimensional "hot spots" 
model 21] for an electron moving in the random field 
of pseudogap fluctuations (considered to be static and 



Gaussian) with dominant scattering momentum trans- 
fers of the order of characteristic vector Q = (7r/a,7r/a) 
(a is the lattice spacing), we use d flH EH the recursion 
procedure proposed in Refs. [22l |23|] . which is controlled 
by two main physical characteristics of the pseudogap 
state: A (pseudogap amplitude), which characterizes the 
energy scale of the pseudogap, and k = £ _1 - the inverse 
correlation length of short range SDW (CDW) fluctua- 
tions. Both parameters A and £, determining pseudogap 
behavior, can in principle be calculated from the micro- 
scopic model at hand [la ]. 

Weakly doped one-band Hubbard model with repul- 
sive Coulomb interaction U on a square lattice with 
nearest and next nearest neighbour hopping was nu- 
merically investigated within this generalized DMFT+E 
self-consistent approach, as described in detail in Refs. 

0E1E1. 

Briefly, the DMFT+E self-consistent loop looks like 
as follows. First we guess some initial local (DMFT) 
electron self-energy E(e). Second we compute the p- 
dependent "external" self-energy E p (e) which is in gen- 
eral case a functional of E(e). Then neglecting interfer- 
ence effects between the self-energies (which in fact is 
the major assumption of our approach) we can set up 
and solve the lattice problem of DMFT [HI El- Fi- 
nally we define effective Anderson single impurity prob- 
lem which is to be solved by any "impurity solver" (we 
mostly use numerical renormalization group - NRG) to 
close DMFT+E equations. 

The additive form of self-energ y is in fact an ad- 
vantage of our approach [3, EE HI- It allows one 
to preserve the set of self-consistent equations of stan- 
dard DMFT [3, El- However there are two distinc- 
tions from conventional DMFT. During each DMFT it- 
eration we recalculate corresponding p-dependent self- 
energy E p (/i, £, [E(w)]) via an approximate scheme, tak- 
ing into account interactions with collective modes or or- 
der parameter fluctuations, and the local Green function 
Gn(iu>) is "dressed" by E p (e) at each step. When input 
and output Green's functions (or self-energies) converge 
to each other (with prescribed accuracy) we consider the 
obtained solution to be selfconsistent. Physically it cor- 
responds to the account of some "external" (e.g. pseudo- 
gap) fluctuations, characterized by an important length 
scale £, into fermionic "bath" surrounding the effective 
Anderson impurity of the usual DMFT. Both cases of 
strongly correlated metals and doped Mott insulators 
were considered [HI, El- Energy dispersions, quasiparti- 
cle damping, spectral functions and ARPES spectra cal- 
culated within DMFT+E, all show a pseudogap effects 
close to the Fermi level of quasiparticle band. 

In Ref. [3] this DMFT+E procedure was generalized 
to calculate two-particle properties, such as dynamic con- 
ductivity, using previously developed recursion procedure 
for vertex corrections due to pseudogap fluctuations [24[ , 
producing typical pseudogap anomalies of optical con- 
ductivity and dependence of these anomalies on correla- 
tion strength U. Below we use the approach of Ref. [II 
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TABLE I: Single-band optical sum rule check in the "hot- 
spots" model. U - dependence. Optical integral in units of 
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FIG. 1: Real part of optical conductivity for strongly corre- 
lated system in the pseudogap state (t' — — OAt, t — 0.25 eV, 
T = 0.0894) in DMFT+E P approximation — U dependence. 
Band filling n = 0.8, pseudogap amplitude A = t, correlation 
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length £ = 10a. Conductivity is given in units of ao = 



to investigate the sum-rule in "hot spots" model. 

To calculate optical integral W we have just used the 
conductivity data of Ref. [19| (extended to a wider fre- 
quency range needed to calculate W), while the r.h.s of 
([3]) was recalculated, using recursion relations for S p (e) 
and the whole self-consistency DMFT+E loop. All calcu- 
lations have been done for a tight-binding "bare" spec- 
trum on the square lattice, with the nearest neighbor 
transfer integral t and next nearest neighbor transfer in- 
tegral t'. 

In Fig. [T]we present our typical data for the real part 
of conductivity (t' = -OAt, t = 0.25 eV, band filling 
n = 0.8, temperature T — 0.089t) for different values 
of Hubbard interaction U = 4t, 6t, lOt, 40t and fixed 
pseudogap amplitude A = t (correlation length £ = 10a). 
It is obvious from these data, that optical integral W is 
different for all of these curves, actually its value drops 
with the growth of U (along with damping of pseudogap 
anomalies [19]). However, the single band optical sum- 
rule §3§ is satisfied within our numerical accuracy, as seen 
from Table I. The small "deficiency" in the values of W in 
Table I is naturally due to a finite frequency integration 
interval over conductivity data of Fig. [T] 

In Fig. [5] we show the real part of optical conductiv- 
ity for doped Mott insulator (fixed U = 40t, t' = -OAt, 
t = 0.25 eV, band filling n = 0.8, T = 0.089£) for dif- 
ferent values of pseudogap amplitude A = 0, A = t, 
A = 2t. Correlation length is again £ = 10a and band 
filling factor n — 0.8. The "violation" of sum-rule here 
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FIG. 2: Real part of optical conductivity for doped Mott 
insulator (U = 40t, t' = ~0At, t = 0.25 eV, T = 0.089t) in 
DMFT+Ep approximation for different values of pseudogap 
amplitude A = 0, A = t, A — 2t. Correlation length £ = 10a, 
band filling factor n — 0.8. 



is especially striking — optical integral obviously drops 
with the growth of A. However, again the single band 
optical sum-rule is strictly valid, as seen from Table 
II. 

To study the details of sum-rule "violation", i.e. the 
dependence of the optical integral W on the parameters 
of the model, we performed extensive calculations of the 
appropriate dependences of the r.h.s. ofEq. ([3]) and opti- 



TABLE II: Single-band optical sum rule check in the "hot 
spots" model. A - dependence. Optical integral in units of 
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FIG. 3: Dependence of normalized optical integral on correla- 
tion strength U in the pseudogap state. All other parameters 
are listed in the figure. At the insert — correlation length 

2 

dependence of optical integral in units of j-t . 



cal integral W on the temperature T, doping, pseudogap 
amplitude A, correlation length of pseudogap fluctua- 
tions £ = k _1 and correlation strength U. Some of the 
results are presented in Figs. [3] -[5] 

Typical dependence of the (normalized) optical inte- 
gral on correlation strength U is shown in Fig. [3] for two 
values of A. We can see rather significant drop of W 
with the growth of U. As to correlation length depen- 
dence, which is shown at the insert Fig. [3j it was found 
to be very weak (practically negligible) in a whole region 
of realistic values of £, so that we shall not discuss it 
further. Dependence of W on the pseudogap amplitude 
A (for several values of U) is shown in Fig. |4j Typical 
doping dependence, which reflects just the dependence 
of the square of the plasma frequency on doping, is 
given in Fig. [5J In all cases other, the change of the rel- 
evant parameters of the model lead to rather significant 
drop in the values of W. As to the temperature depen- 
dence (shown in the insert on Fig. [5J it is rather weak, 
quadratic in T and quite similar to that found in Refs. 



Basically these results show, that the value of the op- 
tical integral depends on all the major parameters of the 
model and, in this sense, its value is not universal, so 
that the optical sum rule is significantly "violated" , if we 
restrict ourselves with a single-band contribution. 



FIG. 4: Dependence of the normalized optical integral on the 
pseudogap amplitude A. The other parameters are listed in 
the figure. 

B. Disordered Anderson — Hubbard model 

In Ref. [2(| we have applied DMFT+S approximation 
to calculate the density of states, optical conductivity and 
phase diagram of strongly correlated and strongly dis- 
ordered paramagnetic Anderson-Hubbard model, with 
Gaussian site disorder. Strong correlations were ac- 
counted by DMFT, while disorder was taken into ac- 
count via the appropria te g eneralization of self-consistent 
theory of localization (2f| [H, We considered 

the three-dimensional system with semi-elliptic density 
of states. Correlated metal, Mott insulator and corre- 
lated Anderson insulator phases were identified via the 
evolution of density of states and dynamic conductivity, 
demonstrating both Mott-Hubbard and Anderson metal- 
insulator transitions and allowing the construction of 
complete zero-temperature phase diagram of Anderson- 
Hubbard model. 

For "external" self-energy entering DMFT+E loop we 
have used the simplest possible approximation (neglect- 
ing "crossing" diagrams for disorder scattering), i.e. just 
the self-consistent Born approximation, which in the case 
of Gaussian site energies disorder takes the usual form: 

£( £ ) = A 2 ]TG( £ ,p) (6) 
p 

where A denotes now the amplitude of site disorder. 

Calculations of optical conductivity are considerably 
simplified [2(| due to the fact, that there are no con- 
tributions to conductivity due to vertex corrections, de- 
termined by local Hubbard interaction. Finally, con- 
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FIG. 5: Dependence of the normalized optical integral on hole 
doping in the pseudogap state. At the insert — temperature 
dependence. All other parameters are listed in the figure 



TABLE III: Single-band optical sum rule check in Anderson 
- Hubbard model. A - dependence. Optical integral in units 
of ^-D. 
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ductivity is essentially determined by the generalized 
diffusion coefficient, which is obtained from the appro- 
priate gene ralization of self-consistent equation of Refs. 
[25l I26U27I |28| . which is to be solved in conjunction with 
DMFT+E loop. 

In Fig. [6]we show typical results for the real part of dy- 
namic conductivity of a correlated metal described by the 
half-filled Anderson-Hubbard model (with bandwidth 
2D) for different degrees of disorder A, and U = 2.5D, 
and demonstrating continuous transition to correlated 
Anderson insulator with the growth of disorder. 

Here again the direct check shows that the single band 
optical sum-rule ([3]) is obeyed within our numerical ac- 
curacy, as seen from Table III. At the same time, optical 
integral W itself obviously changes with disorder. 

Again, to study the details of this sum-rule "violation" , 
i.e. the dependence of W on the parameters of Anderson- 
Hubbard model, we performed detailed calculations of its 
dependences on the temperature T, disorder amplitude 
A and correlation strength U. Some of the results are 



FIG. 6: Real part of dynamic conductivity for half-filled 
Anderson-Hubbard model for different degrees of disorder A, 
and U = 2.5D, typical for correlated metal. Lines 1,2 are 
for metallic phase, line 3 corresponds to the mobility edge 

(Anderson transition) , lines 4,5 correspond to correlated An- 

2 

derson insulator. Conductivity is in units of 

presented in Figs. [7]-[9l 

In Fig. [7]we show the dependence of normalized optical 
integral on Hubbard U , for different degrees of disorder 
(both for strongly disordered metal and correlated An- 
derson insulator). It is seen that in all cases the growth 
of correlation strength leads to rather sharp drop of W 
in metallic state, which becames mush slower in Mott 
insulator. 

In Fig. [8] we present similar dependences on disorder 
strength A. In metallic state optical integral generally 
drops with the growth of disorder, while an opposite be- 
havior is observed if we start from Mott insulator (both 
obtained with the growth of U from metallic state and 
under diminishing U in hysteresis region of the phase dia- 
gram Note the absence of any significant changes in 
the immediate vicinity of critical disorder A C /2D = 0.37, 
corresponding to Anderson metal - insulator transition. 
At the same time it should be noted that the most sig- 
nificant growth of the optical integral takes place as the 
system transforms into disorder induced metallic state, 
obtained from Mott insulator, as observed in Ref . [20I] . 

In Fig. [9] we show the temperature dependence of 
the normalized optical integral, for different degrees of 
disorder. In Anderson - Hubbard model it appears to 
be significantly stronger, than in "hot spots" model (see 
above), and decreases with the growth of disorder. More- 
over, while in relatively weakly correlated state it is qual- 
itatively the same - optical integral diminishes with the 
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FIG. 7: Dependence of the normalized optical integral on 
correlation strength in Anderson-Hubbard model for different 
degrees of disorder A (1,2 - strongly disordered metal, 3 - 
correlated Anderson insulator). 



FIG. 9: Temperature dependence of normalized optical inte- 
gral in Anderson-Hubbard model for different degrees of dis- 
order. At the insert - similar dependence at fixed disorder, 
but for different values of Hubbard interaction U, line 3 here 
corresponds to disordered Mott insulator. 
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FIG. 8: Disorder dependence of the normalized optical inte- 
gral in Anderson-Hubbard model for different values of Hub- 
bard interaction U. Lines 1,2,3 - correlated metal, transform- 
ing into Anderson insulator. Line 4 - Mott insulator state 
obtained with the growth of U from correlated metal, line 5 
- Mott insulator obtained with diminishing U in hysteresis 
region of the phase diagram. 



growth of T, it actually grows in disordered Mott insula- 
tor, as seen from line 3 at the insert in Fig. 9. 

Again, as in the case of the pseudogap "hot spots" 
model, these results for Anderson - Hubbard model 
clearly demonstrate the value of the optical integral is 
not universal and depends on all the major parameters 
of the model and, so the single band optical sum rule is 
strongly "violated". 



III. CONCLUSION 

Based on DMFT+S approach, we have studied the 
single band optical sum rule for two typical strongly 
correlated systems, which are outside the scope of the 
standard DMFT: (i) the "hot spots" model of the pseu- 
dogap state, which takes into account important non- 
local correlations due to AFM(CDW) short-range order 
fluctuations and (ii) Anderson-Hubbard model, which 
includes strong disorder effects, leading to disorder in- 
duced metal-insulator (Anderson) transition, alongside 
with Mott transition. 

We have explicitly demonstrated that the single band 
optical sum rule is satisfied for both models, con- 
firming the self-consistency of DMFT+E approach for 
calculation of two-particle properties. 

However, the optical integral W = 2 J °° Rea(ui)du>, 
entering the single band sum rule ^ is non universal and 
depends on the parameters of the model under consider- 
ation. Most of the previous studies addressed its (rela- 
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tivcly weak) temperature dependence. Here we have ana- 
lyzed dependences on essential parameters of our models, 
showing that these may lead to rather strong "violations" 
of the optical sum rule. As most of the parameters un- 
der discussion may be varied in different kinds of experi- 
ments, these dependences should be taken into account in 
the analysis of optical experiments on strongly correlated 
systems. 
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